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Abstract 

We perform a fuU next-to-leading analysis of the the available experimental 
data on the polarized structure function gi of the nucleon, and give a precise 
determination of its singlet axial charge together with a thorough assessment 
of the theoretical uncertainties. We find that the data are now sufficient to 
separately determine first moments of the polarized quark and gluon distributions, 
and show in particular that the gluon contribution is large and positive. 
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The polarized structure function gi{x, Q^) of the nucleon has been recently measured 
with good accuracy for proton and deuteron targets over a reasonably wide range of 
values of x. This opens up the possibility of a precise determination of the first moments 
of gi, which are of direct physical interest, being related to the nucleon matrix elements 
of axial currents. An extraction of the moments of gi from the data, however, requires 
theoretical input, not only because the data cover a limited range in x but, more impor- 
tantly, because data are obtained at different values of for each x bin, and have thus 
to be evolved to a common value of using the Altarelli-Parisi equations before the 
moments can be evaluated (for recent reviews on the phenomenology of polarized struc- 
ture functions see ref. [§]). We have recently shown Q that these perturbative evolution 
effects can actually be quite large and substantially affect the extraction of the moments 
from gi. Furthermore, effects which are formally next-to-leading order (NLO) may lead to 
significant evolution because of the large contribution of the polarized gluons to gi driven 
by the axial anomaly. 

The recent computation of the full matrix of two-loop anomalous dimensions makes 
a consistent NLO analysis of gi now possible. It is the purpose of this paper to perform 
such an analysis, and use it to provide a precise determination of the first moment of 
gi. After reviewing the NLO formalism, discussing scheme dependence, and the effect of 
NLO corrections on the LO small-x behaviour of parton distributions, we will use it to 
extract polarized parton distributions from the data. We will find that the data allow 
a determination of both the quark and gluon distributions without the need for extra 
theoretical assumptions, and in particular strongly constrain their overall normalizations 
and small-a; behaviours: we will thus be able to infer the existence of polarized gluons in 
the nucleon from an analysis of the scale dependence of gi. We will then use these parton 
distributions to determine the first moment of the structure function gi and the nucleon 
matrix element of the singlet axial current, or singlet axial charge, whose unexpected 
smallness has attracted a good deal of interest. We will finally provide an evaluation of 
the various sources of statistical and systematic error related to these determinations. 

The structure function gi is related to the polarized quark and gluon distributions by 

gi{x, Q2) = [Cm ® Aqns + Cs ® AS + 2nfCg ® A^] , (1) 

where (e^) = ^ X]r=i ' ® denotes the usual convolution with respect to x, the nonsinglet 
and singlet quark distributions are defined as 

A?NS = ( (5) - ^^^^ ^ " 5Z(A?^ + A^z), (2) 

1=1 i=l 

where Aqi and Aqi are the polarized quark and antiquark distributions of fiavor i, and 
Ag is the polarized gluon distribution (see ref. Q for further details of notations and 
conventions.). The polarized parton distributions evolve according to the Altarelli-Parisi 
equations 

at iTx , . 

rf/AE\_a.(t)/ Pi, 2n,ps\ /AE^ 
dt\Ag)- 2. \ )®\Ag 
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where t = ln((5^/A^). The coefficient functions C and sphtting functions P may each be 
expanded in powers of dg: at NLO 



C{x,as) = C^^Hx) + —C^^Hx) + 0{al) (4) 
P{x, as) = P^'\x) + ^P^'Hx) + 0{al). (5) 

In accordance with the partonic picture C^°g(a;) = Cq^\x) — S{1 — x), while Cg^\x) = 
0. It will also prove convenient to introduce anomalous dimensions 7(A^, cts) = 
Jq (ix x^~^P(x, cts), i.e. the Mellin transforms of the splitting functions, as well as anal- 
ogously defined moment-space coefficient functions C{N,as) and parton distributions 
AQNs(iV,Q'), AE(iV,Q2) and Ag{N,Q^). 

Whereas the NLO coefficient functions C^^^ have been known for some time [0 the 
two loop splitting functions P^^^ have been only recently determined 0,0 only their first 
moments being known previously 0. The NLO coefficient functions may be modified 
by a change of factorization scheme which is partially compensated by a corresponding 
change in the NLO anomalous dimensions, hence both are required in order to specify a 
NLO computation completely. Previous analyses which included NLO effects only in the 
coefficient functions (such as [Q) were thus necessarily incomplete and treated only the 
first moments consistently at NLO. It is now possible to test explicitly whether, as was 
claimed in ref. 0], the NLO gluon contribution to the first moment of gi is the dominant 
NLO effect, and furthermore whether it has a sizable effect on the dependence of gi. 

The NLO anomalous dimensions and coefficient functions of ref. are given in the MS 
scheme. Since chiral symmetry is respected, matrix elements of nonsinglet axial currents 
are conserved, nonsinglet axial charges do not evolve, and thus AgNs(l, Q^) is independent 
of Q^. In all such schemes at NLO 

Cns(1, as) = Cs(l, as) = l- — + 0{al) (6) 

n 

(although at higher orders CNs(l,as) 7^ Cs{l:0:s))- Matrix elements of the axial singlet 
current are instead not conserved because of the axial anomaly, so that the singlet axial 
charge depends on scale. In the MS scheme the first moment of Cg^^ vanishes, the gluon 
decouples from the first moment of gi and the scale dependent singlet axial charge is thus 
equal to AJ:{1,Q'^). 

Factorization schemes where this happens are somewhat pathological, in that they 
include soft contributions to the cross section in the coefficient function rather than ab- 
sorbing them completely into the parton distributions []8|-p!0| . We will instead perform our 
calculations in schemes where all soft contributions are properly factorized into the parton 



distributions. The first moment of the gluon coefficient function at NLO is then [jTT 

Cg{l,as) = -^ + Oial). (7) 



With our conventions the expressions of 7^°^ and 7^^^ of ref. must be divided by —4 and 
-8 respectively and 'jqg should be further divided by 2n/. 



2 



With this choice, the first moment of gi is given at NLO by 



ri(Q2) ^ [\,(x,Q')dx = ^ [(1 - ^){Aq^s{l.Q^) + m^.Q^)) -nf^Ag{l,Q^)] , 
Jo 

(8) 

and the first moment of the scale dependent eigenvector of the singlet Altarelli-Parisi 
equations (^) is 



aoiQ') = AE(1, Q2) _ Uf^Agil, Q^). (9) 

Ztx 

The corresponding eigenvalue of the anomalous dimension matrix coincides with the 
anomalous dimension of the singlet axial current, so ao is identified with the singlet axial 
charge 

{p,s\j^\p,s) = Ms^'ao{Q^), (10) 

where p, M and s are the momentum, mass and spin of the nucleon. The other eigenvector 
of perturbative evolution is the first moment of the polarized singlet quark distribution, 
AE(1,(5^), which is then independent of , and may be identified with the conserved 
singlet quark helicity |Tl| - p!3[] . 

In fact the eigenvectors remain the same to all orders in perturbation theory, because 
of the Adler-Bardeen theorem [^, which states that the NLO mixing of the divergence of 
the singlet axial current with a gluonic operator (the anomaly), which is responsible for 
its scale dependence, does not receive higher order corrections. Thus if we require that 
AE(1, Q"^) is scale independent then the axial charge is given by eq. (^ to all orders. This 
means that the first moments of the singlet quark and gluon coefficient functions are not 
actually independent: to all orders in perturbation theory (||) becomes simply 

ri(Q2) = [CNs(l,«.)AgNs(l,Q') +C^s(l,«s)ao(Q')] , (11) 

and thus 

(y 

C,(l,a,) = -^Cs(l,a.). (12) 

Given anomalous dimensions and coefficient functions in a particular factorization 
scheme any other factorization scheme can be constructed ||15|| by introducing a scheme 
change specified by a function Z]<iQ{N,as) = 1 + f^^Ns(-^) ^("^s) a matrix 

ZQ,{N,aa) = 1 + ^z^\n) + 0{a1). The NLO anomalous dimensions and coefficient 
functions then change according to 

7^.sHAr)-.7(i)(iV)-f.W(iV), 

- l^i\N) + [z^i\N),^f{N)] - ^4\n), (13) 

C^s\n)^C^^\n)-z^I){N), (14) 
where (3q = 11— "^Uf is the one loop coefficient of the QCD beta function. 
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For simplicity we will only discuss scheme changes where zns = Zqq, i.e. such that the 
relative normalization of the singlet and nonsinglet quark distributions is unaffected. The 
conservation of the nonsinglet axial current then fixes the first moment ^^^g = z^q — 
whenever the original scheme respects chiral symmetry. We then wish to consider specif- 
ically factorization schemes in which the first moments of the coefficient functions satisfy 
eq. (1^): starting from the MS scheme, eq. (|4|) then fixes z^g {1) = 1.11 In order to com- 
pletely specify the first moment of z we use the Adler-Bardeen condition that the two-loop 
eigenvector of perturbative evolution as given by (||) be identified to all perturbative orders 
with the matrix element of the axial current. Knowledge of the NNLO anomalous dimen- 



sion of the axial current then fixes the first moments of the remaining two entries of 
the matrix z. 

We will consider several schemes which differ in the way the remaining moments of 
the coefficient functions and anomalous dimensions are constructed. In the first scheme, 
we simply take z{x), the inverse Mellin transform of z{N), to be independent of x. This 
scheme is thus the minimal modification of the MS scheme such that the first moments 
of parton distributions satisfy the anomaly constraint eq. (j^); in particular, the large and 
small X behaviour of the coefficient functions and anomalous dimensions are then the same 
as in MS. We will refer to this as the Adler-Bardeen (AB) scheme. The matrix which 
transforms from MS to the AB scheme is 

where, for SU(3) color, Cf = |, Ca. = 3 and Tp = ^. Notice that the two lower entries in 
the MS NLO anomalous dimension matrix ^ turn out to be already consistent with the 
Adler-Bardeen condition above, and NNLO anomalous dimensions Jl^; the corresponding 



entries of the scheme change matrix eq. (|T^) therefore vanish. 



Transformations such as (jT^) which take us from MS to a scheme where the gluon 
contributes to the first moment of gi correspond to removing soft contributions from the 
coefficient functions |l9|-|T0[]. Rather than doing this by hand, as in the AB scheme above, 
the subtraction may be performed by computing the coefficient functions in the presence 
of an explicit infrared regulator, which automatically enforces eq. Jl^. The entries Zqq 



and Zqg of the z matrix are then fixed using eq. (|T^) and the MS coefficient functions 



V)|ms = (i - NW+T) + Si{N)) - S2{N) 



2N'^{N+1) 



(16) 



where Sj{N) = X]fc=i W- '^^^ lower entries of the transformation matrix can then be 
taken to be zero as in the AB scheme. 



^ A scheme change of this kind was constructed in ref. [16|; the form of the matrix z given 
there appears however to be incorrect. Also, note that the partial result for the NLO splitting 
functions given there is incorrect, as explained in ref. 
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One possibility is to renormalize while keeping the incoming particle off-shell (OS 
scheme, henceforth); the coefficient functions are then given by 



^9 )\os ~ N^N+iy-' ■ 



(17) 



An alternative option is to endow the quarks with a finite mass (Altarelli-Ross, or AR 
scheme); the quark and gluon |11,10| coefficient functions are then 



r^s(l)^Anl _ rp N'' + l+N(N-l)Si(N) 
^9 >\aR - N-^N+1) ■ 

(18) 

Notice that in all three of these schemes the first moments of the coefficient functions are 
given by eqs. (^ and (|^), and thus the NLO relation between the first moment of gi and 
the singlet axial charge ao implicit in eqs. (H) and (||) is automatically satisfied. 

The main effect of the NLO corrections to perturbative evolution is the coupling of the 
gluon to and in particular its contribution to the first moment Fi eq. (||), which does 
not decouple as oo [|ll|. However, NLO corrections may also substantially affect 



the small- a; behaviour of parton distributions and coefficient functions. Indeed, unlike the 
unpolarized case, the NLO contributions to the polarized splitting functions and coefficient 
functions display a stronger singularity as x — > than their LO counterparts; accordingly 
their Mellin transforms display a stronger singularity as ^ 0. More specifically the 
singularities in the NLO MS anomalous dimensions and coefficient functions take the 
form 

^irsV) = ^ {'^CaCf - 3Cl) + O , (19) 

(^)(m = — ( -^CFTFnf + 2CACF-^Cl -2CaTf-CfTf \ ±_ 
A3 V 4CaCf + 2C2 8Cl~4CFTFnf J^^ym 

- Ci^^ + --Tf^ + O (1) , (20) 

whereas the LO anomalous dimensions have a simple pole in A (details of which may be 
found in ref. Wi)E 



3 Because the quark mass breaks chiral symmetry, it is now necessary to perform an extra 
subtraction in order to ensure that the nonsinglet axial currents are conserved. 

^ The presence of double logarithms in the NLO polarized splitting functions and coefficient 
functions strongly suggests that the systematic cancellation of collinear singularities which char- 
acterizes the small-x behaviour of unpolarized splitting functions does not occur in the polarized 
case. In the nonsinglet channel a summation of these double logarithmic singularities to all orders 



in as has been attempted in ref. |1S]. 
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The NLO corrections could therefore have a significant impact in the smaU x region, 
and in particular require a summation of logarithmic effects in ^, which could be done in 
analogy to the unpolarized case ||2^ if the coefficient of the most singular contributions to 
the polarized splitting functions as x — > were known to all orders in as- In order to assess 
at which values of x and these effects might begin to be relevant (and in particular 
whether they already amount to a sizable correction in the presently measured region) gi 
may be determined by solving the NLO evolution equations with the approximate small- 
form eq. (|19D of the anomalous dimensions: 



AgNs(iV,Q') =A?Ns(iV,Qo) 



v±{N,Q^) =v±{N,Ql) 





^0 







1 + 



iV3 



(21) 



1 + 



iV3 



where v± and \± are the eigenvectors and eigenvalues of the LO singlet anomalous dimen- 
sion matrix 7*^*^) (see ref. [§[), as{Q'^) is computed at NLO, Qq is the starting scale, and 
the coefficients e are explicitly given by 



37I-/30 ' 



112 

37r/3o 



(l-^)±i(l 



1 _ '^HL 
^ 32 



(22) 



The NLO corrections do not mix the LO small x eigenvectors, because mixing terms are 
O (772)- The eigenvectors of perturbative evolution at NLO are thus the same as at LO: 
at small x and large Ag and AE have opposite sign 0, and in particular (for any 
plausible parton distributions) AE < and Ag > 0. It is interesting to observe that this 
result is scheme independent, because the leading singularities in the NLO corrections 
only receive contributions from the diagonal projections of the NLO anomalous dimension 
matrix onto the LO eigenvectors, which are themselves scheme independent. 

The leading NLO small-x behaviour can now be found by inverse Mellin transform of 
eq. (^1]) in the saddle point approximation: 



AgNs(a;,Q^ 



l + 6Ns(r^)'(«.(Q^)-«.(Q^)) 



7ns 
3 



(23) 



1 a^s(Ql) ■ 



where J\f are normalization constants, a = p = a/^/C, ^ = C 
a reference value of x such that the approximate small-x form of the anomalous dimensions 
is applicable for x ^ xq and ^ Qq, and 7^g = while 7± are as given in ref. 



^ Eq. ( p3| ) only gives the NLO generalization of the asymptotic LO behaviour a^^^^e^^'^ when 
the boundary conditions are soft, as discussed in ref. If the boundary condition is hard, for 
example e^"^, then for p ^ 7/A the LO behaviour reproduces the boundary condition Q, and 
the NLO correction to it is given by + tiX^^ (qs((5o) ~ 0(s{Q^))~\ , where ti = cns, e± in in the 
nonsinglet and singlet cases respectively. In this case the NLO correction is thus x-independent. 
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The terms in square brackets in eq. (|23| ) give the NLO correction to the LO asymptotic 
small X behaviour P,|2T[] . Because all coefficients e eq. ([2^ ) are positive the NLO corrections 



lead to a further increase proportional to ^^/^ of the parton distributions at small x. The 
coefficient of this increase is however rather small, for instance with Uf — A one gets 
cns/Tns ~ e+/7+ ~ so that the correction is small in the presently accessible small x 
region. These conclusions however only apply to the region where the no summation of 
logs of ^ is necessary so that the NLO in as may be treated as a subleading correction, 
and could be substantially altered at smaller values of x. 

The leading small x behaviour of gi can be found at NLO using the small x NLO 
solution eq. (|2T| ) and coefficient functions eq. (|20| ) in the Mellin transform of the expression 
of gi eq. (|l|). Because the NLO coefficient functions (|20|) only have a singularity at NLO 
they actually do not contribute to the leading small-x behaviour of gi, which is thus found 
by simply taking the appropriate linear combination of the small- a; parton distributions 
eq. (|2^). The NLO correction to the small-x behaviour of the coefficient functions may 
nevertheless have a significant impact, especially at low scales (i.e., when Q is close to 
the starting scale Qq) where evolution effects are negligible. Indeed, the singularity 
corresponds to a log - rise of the coefficient function, and would therefore lead to a rise of 
both singlet and nonsinglet contributions to the structure function gi even if the parton 
distributions themselves did not rise. The coefficient of this rise is however not scheme 
independent: for instance, the coefficient of the leading singularity in the gluon coefficient 
function is the same in the MS, AB and AR schemes, but is twice as large in the OS 
scheme. 

Having established that a NLO treatment is adequate in the region of current exper- 
imental data, we can proceed to a determination of the physical observables related to gi. 
Even though our purpose here is not to establish a parametrization of polarized parton 
distributions, we have to construct such a parametrization since only LO parametrizations 
are currently available^. We parametrize the initial parton distributions according to 

A/(x, Ql) =N{af, /3f,af) ry^x"^ (1 - xff (1 + a/x), (24) 

where N{a,f3,a) is fixed by the normalization condition N{a,P,a) j}dxx'^{l — x)^(l -|- 
ax) = 1, and A/ denotes AE, Aqns or Ag. In our previous analysis g], which used only 
proton data, the respective three sets of parameters could not be independently determined, 
while the small x behaviour had to be fixed and then varied in a plausible range. Now, by 
including the deuteron data 0, we can disentangle the nonsinglet and singlet quark and 
gluon contributions to gi, because gi is dominated by the isotriplet term, which contributes 
about 90% of its first moment, while gf is isosinglet. We can thus determine independently 
almost all of the parameters of the singlet and nonsinglet parton distributions (including 
the parameters o;/, which determine their small-x behaviour). 

We determine gi at all x and by solution of the NLO evolution equations with 
boundary conditions of the form ( |2^ ) at Qq = 1 GeV^. The various parameters in eq. 



^ A comprehensive review of the present status of polarized parton parametrizations is given 
in ref. [ p2[ . A NLO parametrization has been presented in ref. [ p3[ , but in the unsubtracted MS 
scheme, which, as previously discussed, is not properly factorized. 
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are then found by fitting gi{x,Q^) to the recent precision experimental determination of 
gi for proton and deuteron which are given along a curve ~ Q'^{x) for each exper- 
iment. As in gi is extracted from the measured asymmetry Ai using the leading- twist 
expression^ gi{x, Q"^) = Ai{x, Q^); Fi is in turn computed from the SLAG deter- 



mination of the ratio Q ) of the longitudinal to transverse virtual photoabsorption 



cross section, and the most recent NMC determination [25| of F2{x,Q'^). The deuteron 



structure function is defined as the average of the proton and neutron structure functions 
and is obtained from the deuteron asymmetry after applying a correction to account for 



(i-wave admixture 26 



gi{x, Q2) ^ 1 (^P(^, q2) + gn^^^ g2)) ^ ^d^^^ Q^)Ft{x, _ 1.5^^) (25) 

where ood = 0.05. 

The normalization of the nonsinglet quark distribution at Qq (which lies below the 
charm threshold) is fixed by assuming SU(3) symmetry of the matrix elements of the axial 
current, determined from hyperon j3 decays: 



77ns = / AQNs(a;, Q'^)dx = ±\gA + Jas, (26) 
= 1.2573 ±0.0028; ag = 0.579 ± 0.025, (27) 

where the plus (minus) sign refers to a proton (neutron) target. While the impact of 
possible SU(3) violation on our results will be discussed below, we will defer to a subsequent 
publication the possibility of using the data themselves to test SU(3) or SU(2) (and thus the 
Bjorken sum rule) . At higher scales, new nonsinglet contributions arise as the various heavy 
quark thresholds are passed, so that Aq'Ns(l5 Q^) is not scale independent across thresholds. 
The corresponding heavy quark distributions are generated dynamically, assuming that 
they vanish on threshold, and imposing continuity of ao(Q^) eq. (||). 

It turns out that the data are good enough to determine all the remaining parameters 
directly with a single exception, namely, the parameters and which control the shape 
of the singlet distributions AE and Ar? at intermediate x, and which are difficult to pin 
down individually since the evolution mixes AE and A^f. We will thus take = a^, 
and determine the remaining ten parameters from a fit to the data. The respective best 
fit values are listed, in table 1 for the AB, AR and OS schemes. The errors given there 
are statistical errors from the fit, computed by taking into account correlations. The 
corresponding determination of gi{x) in the AB scheme is shown in fig. 1 at the starting 
scale as well as at the scale Q{x) of the various data sets and at = 10 GeV^. 

The main features of this determination of polarized parton distributions are the fol- 
lowing. The evolution of gi is rather similar to that discussed in ref. , thus demonstrating 
that the main NLO effect is indeed the direct gluon contribution to the first moment of gi 
which was already included there; the main difference between the present determination 



^ Kinematic higher twist corrections which are sometimes included are neglected here for con- 
sistency since no systematic treatment of these corrections is available. 
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of polarized parton distributions and that of ref. [|| is the more detailed nature of the fit, 
due to the inclusion of the deuteron data. 

Specifically, the large x behaviour, controlled by the exponent /3, is in fair agreement 
with the expectations based on QCD counting rules ||2^, which predict (3q — /3ns — 3 and 



Pg ~ 4; of course these parameters are strongly scheme dependent (as explicitly displayed 
in table 1). The value of (3g cannot be determined very accurately since large-x data are 
at relatively large Q^, where the direct coupling of the gluon to gi is suppressed by the 
small value of as] in the OS scheme the results are largely independent of this coefficient 
which thus cannot be fitted at all. 

The small x behaviour turns out to be at least qualitatively constrained by the data. 
The initial singlet parton distributions are found to display soft valence-like behaviour, i.e. 
not to grow at small x, as was predicted by Regge theory [^^. However the nonsinglet 



quark distribution is found to grow approximately as Aq ~ -j=B We have explicitly 

verified that these behaviours indeed correspond to an overall global minimum, and that 
the quality of the fit deteriorates uniformly as the values of a/ deviate from those given in 
table 1, by repeating the fit with fixed values of the three exponents ct/ chosen in the range 
—0.9 < CK/ < +0.5. The singular behaviour of the nonsinglet quark is the main difference 
between these results and those discussed in ref. where, since it was impossible to 
determine it from the data, it was assumed to be valence-like. 

Even though these results are suggestive, they should be taken with some care: firstly, 
the values of a/ are strongly scheme-dependent, and partly just refiect the shape of the 
coefficient functions at small x and small in the various schemes: for instance, the fact 
that the coefficient of the leading singularity in Cg^"^ is larger in the OS scheme explains 
why the value of ag in this scheme is larger. Furthermore, these are not necessarily the 
asymptotic behaviours of the various parton distributions as a; 0: in particular, the 
large value and uncertainty in ons indicate that the asymptotic behaviour of gf is still 
setting in at the smallest experimental x values; likewise, gf at small and small x grows 
large and positive with the present parametrization, but this happens outside the data 
range (the growth is barely seen starting on the = 1 GeV^ curve of Fig. lb). In this 



respect, present-day polarized data are to be compared to the NMC unpolarized data p5[] : 
a kinematic coverage comparable to that available at HERA for determining the small- 
X behaviour F2 would be required to determine precisely the small x behaviour of polarized 
parton distributions. This is reffected by the large uncertainties in the parameters a/ which 
govern the transition to the small x region. 

Finally, the data turn out to allow a good determination of the size of the quark 
and gluon distributions. The size of the gluon distribution drives perturbative evolution; 
the evolution already observed due to the fact that the SMC and E143 data are taken at 
different values of for equal x (displayed in fig. 1) turns out to be sufficient to require rjg 
(which gives the first moment of Ag at the initial scale Q = Qq) to be large and positive. 
The relation eq. (^) between the first moments of the quark and gluon distributions and the 
first moment of gi then inevitably leads to a rather large quark singlet. This is a remarkable 



This seems to be qualitatively consistent with the summation of double logarithms presented 
in refs.| p^j30| ], since this suggests that both the polarized and unpolarized nonsinglet quark 
distributions have the same leading small x behaviour, at least in perturbation theory. 
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result: the scale-independent first moment of AE, given by 77^, is found empirically to be 
equal to ag eq. within errors. The large violation of the Zweig rule in the first moment 
of gi appears then to be almost entirely due to a large perturbative gluon contribution to 



gi, as was conjectured in ref. |11 



A large gluon distribution implies substantial evolution effects, and thus a substantial 
correction due to the determination of the moments of gi from the experimental data 0] . 
Indeed, we can now determine the first moment eq. (|^), as well as the singlet axial charge 
of the nucleon defined according to eq. (]^) from the best-fit polarized parton distributions. 
The results in the AB scheme, which we shall take as a baseline, are displayed in the first 
row of table 2. Only is shown; the deuteron is obtained from it by subtracting the 
isotriplet contribution r(^^ = j^CnsII, tts)fi'A since we have assumed isospin to be exact. 
The value of Fi is found to be lower than the value extracted neglecting the evolution 
of the asymmetry, and consequently the value of the axial charge ao we obtain is about 
a half of the value quoted by the experimental collaborations. It is thus important to 
establish whether this result is a robust consequence of the inclusion of NLO evolution 
effects, or whether it is a by-product of the procedure used so far. Furthermore, if the 
result is confirmed, it inevitably implies an increase in the theoretical uncertainty in the 
determination of the axial charge, which must be more accurately estimated. 

First, we test the sensitivity to the specific functional form eq. (^ip of the parton 
distributions. To this purpose, we change the parameterization (^4]) by replacing the last 
factor with (1 + hf^fx -f- ajx), and then repeating the fit with various choices for the new 
parameters 6/. In particular, we fix hq = hg and then we vary —5 < 6^, 6ns ^ 5. No 
significant variation in the results or the quality of the global fit is found, indicating that 
the fit is stable. The results for extreme values of 6ns while bq = bg = are shown in 
table 2; variations of bq, bg produce even smaller effects. 

Other sources of uncertainty related to the fitting procedure are due to the fact that the 
structure function is not measured over the full range of values of x, and is only sampled 
at a finite number of points. The former uncertainty comes from the extrapolations to 
small and large x, and is thus already included through the propagated errors in the six 
parameters a/ and It is interesting to note that if we were to assume that the initial 
nonsinglet quark distribution was fiat at small x, i.e. q;ns — 0, Fi would scarcely be altered, 
even though rjq and rjg both fall substantially (see table 2). The latter effect, however is due 
more to the deterioration in the quality of the fit at intermediate x than to a large change 
in the small x contribution. The other uncertainty is related to the fact that we determine 
the first moments of gi and the parton distributions by integrating the respective best-fit 
forms. This is to be contrasted with the procedure followed by experimental collaborations, 
which instead determine the first moment of gi by summing the data over the experimental 
bins; evolution effects could then be included as corrections to each bin separately. Such 
a procedure is however problematic because evolution effects are actually rather large in 
many of the bins, and furthermore the value of gi{xo, Qq) is determined (nonlinearly) by 
the measurements at all x > xq and < Qq so it is not possible to disentangle truly 
independent corrections to individual data points. Thus a simple sum over bins ignores 
the constraints imposed by perturbative evolution on the shape of gi in the (x, Q^) plane 
and gives undue weight to particular data points. The integration of fitted distributions 
has nevertheless an inevitable sampling uncertainty. To estimate this we have computed 
the difference between the experimental value of gi in each bin and our best-fit gi{xb, Qf), 
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where (x^, Ql) are the central values of x and for each x bin, and have then computed 
the first moment A of this difference over the measured region by multiplying it by the 
bin width and summing over bins. We find Ap(SMC) = -0.0032, Ap(E143) = 0.0007, 
A'^(SMC) = -0.0045, A'^(E143) = 0.0010 for the four experiments, showing that this 
correction is indeed small and does not lead to a systematic bias. Combining these results 
we obtain an estimate for the overall statistical uncertainty due to the choice of the fitting 
procedure, to be added to that of the fit itself (see table 3). 

We now turn to the various sources of systematic theoretical uncertainty. All the cor- 
responding results, along with the of the various fits, are summarized in table 2. First, 
we consider the possibility of SU(2) or SU(3) violation. The former will in general induce 
current mixing effects: for instance, in the presence of isospin violation, the proton and 
deuteron matrix element of the isosinglet current may be unequal to each other. However, 
due to the expected smallness of isospin violation, this effect is negligible compared to 
the error on the values of qa and ag induced by SU(2) and SU(3) violation in the quark 
distributions, which are then the dominant source of uncertainty. We estimate these by 
assuming a 2% uncertainty on qa (of the same order of the accuracy to which isospin 
symmetry of unpolarized quark distributions may be established [^) and a 30% uncer- 
tainty on ag [^. Notice that when these parameters are varied the size of the nonsinglet 
first moment of gi varies; however, table 2 shows that the best-fit value of the singlet first 
moment and thus of ao adjusts itself and remains surprisingly stable. 

The uncertainty related to the value of as, which is not negligible in a NLO com- 
putation, is simply estimated by repeating the fit when as is varied in the range 
as{Mz) = 0.117 ± 0.005. The position of the quark thresholds is varied from 0.75mg to 
2.5mg (with rric = 1.5 GeV, mf, = 5 GeV). 

An important source of theoretical uncertainty is due to the lack of knowledge of higher 
order corrections, as reflected by the dependence of the results on factorization scale 
and renormalization scale /i^. We estimate this by taking = kiQ'^ and /x^ = k2Q^, 
and varying 0.5 < ki, k2 < 2. When the scales are varied the physical parameters turn out 
to have a stationary point within this range; the associated error is thus asymmetric, but 
rather large, consistent with the fact that evolution effects are important. The fluctuations 
found when varying the factorization scale are consistent with the spread of the results 
displayed in table 1 in various factorization schemes. We have also checked that including 
the known two and three loop corrections to the flrst moments of the coefficient functions 
produces similar variations. 

The errors corresponding to all these sources, given by the maximal variation of the 
results as the parameters are varied in the respective ranges, are summarized in table 3. 
We have not included an error due to higher twist terms since we know of no reliable way of 
estimating it: experience with unpolarized data suggests that it is probably rather smaller 
than the error from higher order corrections. 

In conclusion, the analysis of various sources of theoretical uncertainty confirms the 
stability of our determination of the first moment of gi and the axial charge of the nucleon. 
Adding in quadrature the various sources of error, summarized in table 3, we get finally 

r?(3 GeV^) = 0.118 ± 0.013 (exp.) (th.), 

rf(3 GeV^) = 0.024 ± 0.013 (exp.) +J;JJJ (th.), 

r^(10 GeV^) = 0.122 ± 0.013 (exp.) (th.), (28) 

rf (10 GeV^) = 0.025 ± 0.013 (exp.) +J;JJJ (th.), 

ao(10 GeV^) = 0.14 ± 0.10 (exp.) (th.). 
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The deuterium values refer to the structure function eq. (^); they must be multiphed by 
l — l.bujD in order to compare with the results quoted by the experimental collaborations. 
Values of ao at any other scale are obtained using NLO evolution 



«o(Q')= [l + %^^ + OK)]ao(oo), (29) 

which the expression (^ ao satisfies by construction. The dominant source of theoretical 
uncertainty is that related to higher order corrections, i.e. to renormalization and fac- 
torization scale. The experimental uncertainty includes the various errors related to the 
fitting procedure (hence also the experimental systematics) , and is still dominant. This 
is mainly due to the fact that the data allow us to only partly constrain the shape of the 
polarized gluon distribution, and thus the perturbative evolution of gi. 

It is interesting to compare the value of the axial charge eq. (|28|), obtained using 
eq. (^, with that obtained using eq. (JTj) with the coefficient functions expanded to NLO 
as in eq. (P): 

a'oiQ') ^ [Cs(l,a.)]-^ {j^^iiQ^) - Cns(1, «.)Aqns(1, Q')) • (30) 

If we could work to all orders ao = ag; however, in a k-th order perturbative computation 
the two determinations differ by {k + l)-th order corrections. Indeed, at NLO 

a'oiQ') - aoiQ') = - 2n/ (^)' A^7(l, Q^) + 0(af ). (31) 

This difference may be quite large if the gluon distribution is large: indeed, using eq. (|5^) 
we get ao(10 GeV^) = 0.09 ±0.15 (stat.), which differs considerably from the central value 
eq. (|28|), though within the theoretical error. Notice that ( |3T| ) is scale dependent: while 
it vanishes asymptotically, at low scales it is large enough that ao((5^) actually increases 
with below 10 GeV^ rather than decreasing as the axial charge should, according to 
eq. (^. 

Finally, we can estimate the size of the quark and gluon distributions in the AB 
scheme: ^ 

/ dxAJ:{x) = 0.5 ±0.1, 

■^•^ (32) 
dxAg{x,lGeV^) = 1.5 ± 0.8, 

where all errors have been added in quadrature. The gluon distribution is rather large even 
at the quoted low scale (the corresponding value at 10 GeV^ is roughly twice as large) and 
even though the error is large, it differs significantly from zero. The value of the (scale 
independent) singlet quark distribution is in agreement within errors with the prediction 
of the Zweig rule, which would identify it with ag eq. (^7|) . 

In summary, we have given a NLO determination of the main physical observables 
related to the polarized structure function gi. Our main result is that the data already 
constrain the size of polarized parton distributions, and in particular require a rather large 
polarized gluon distribution. This in turn implies that perturbative evolution effects are 
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not negligible, and in fact substantially affect the extraction of the first moment of gi 
from experimental data. Indeed, we find a value of the nonsinglet first moment which is 
significantly smaller than that obtained from a purely LO analysis. More importantly, 
we show that the error on the determination of the singlet axial charge of the nucleon 
is significantly larger than usually recognized, due essentially to the unknown effects of 
higher order perturbative corrections. The recent precise data on gi are providing us with 
surprisingly accurate information on polarized parton distributions, but they also show 
that the theoretical and phenomenological interpretation of these data is significantly more 
subtle than previously expected. 

Acknowledgement: Wc thank R. Mertig and W. van Neerven for providing details 
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about the experimental data, and G. Altarelli, R.K. Ellis and G. Mallot for discussions. 



13 



References 



[1] E143 Collaboration, K. Abe et al., Phys. Rev. Lett. 74 (1995) 346; 

SMC Collaboration, D. Adams et al., Phys. Lett. B329 (1994) 399. 
[2] E143 Collaboration, K. Abe et al., Phys. Rev. Lett. 75 (1995) 25; 

SMC Collaboration, D. Adams et al, Phys. Lett. B357 (1995) 248. 
[3] G. Altarelli and G. Ridolfi, Nucl. Phys. B (Proc. Suppl.) 39B (1995) 106; 

S. Forte, [hep-ph/9409416| , in "Radiative Corrections: Status and Outlook", 

B. F. L. Ward, ed. (World Scientific, Singapore, 1995); 

R. D. Ball, in the proceedings of the "International School on Nucleon Spin Structure" , 
Erice, August 1995 (to be published). 
[4] R. D. BaU, S. Forte and G. Ridolfi, Nucl. Phys. B 444 (1995) 287. 
[5] R. Mertig and W. L. van Neerven, Leiden preprint INLO-PUB-6/95; Erratum, private 

communication (November 1995). 
[6] G. Altarelh and G. Parisi, Nucl. Phys. B126 (1977) 298. 
[7] J. Kodaira et al., Phys. Rev. D20 (1979) 627; Nucl. Phys. B159 (1979) 99; 

J. Kodaira, Nucl. Phys. B165 (1979) 129. 
[8] R. D. Carlitz, J. C. Collins and A. H. Mueller, P/iys. Lett. B214 (1988) 229. 
[9] G. Altarelh and B. Lampe, Zeit. Phys. C47 (1990) 315. 
[10] W. Vogelsang, Zeit. Phys. C50 (1991) 275. 
[11] G. Altarelh and G. G. Ross, Phys. Lett. B212 (1988) 391. 
[12] S. Forte, Nucl. Phys. B331 (1990) 1. 

[13] G.M. Shore and G. Veneziano, Phys. Lett. B244 (1990) 75; Nucl. Phys. B381 (1992) 
23. 

[14] S. Adler and W. Bardeen, Phys. Rev. 182 (1969) 1517. 
[15] W. Furmanski and R. Petronzio, Zeit. Phys. Cll (1982) 293. 
[16] E.B. Zijlstra and W.L. van Neerven, Nucl. Phys. B417 (1994) 61. 
[17] S. A. Larin, Phys. Lett. B334 (1994) 192. 

[18] M. Stratmann, A. Weber and W. Vogelsang, Dortmund preprint DO-TH 95/15, [hep- 
I ph/9509236 

[19] J. Bartels, B. I. Ermolaev and M. G. Ryshkin, DESY preprint 95-124. 
[20] R. D. Bah and S. Forte, Phys. Lett. B351 (1995) 313. 
[21] M. A. Ahmed and G. G. Ross, Phys. Lett. B56 (1975) 385; 

M. B. Einhorn and J. Soffer, Nucl. Phys. B74 (1986) 714 ; 

A. Berera, Phys. Lett. B293 (1992) 445. 
[22] T. Gehrmann and W. J. Stirling, Durham preprint DTP/95/78, |liep-ph/9510243 . 
[23] M. Gliick et al, Dortmund preprint DO-TH 95-13 |hep-ph/9508347| . 
[24] L. W. Whitlow et al., Phys. Lett. B250 (1990) 193. 



14 



[25] NMC Collaboration, M. Arneodo et al. preprint CERN-PPE/95-138, [Eep^ 

I ph/9509406 . 

[26] L. L. Frankfurt and M. Strikman, Nucl. Phys. A405 (1983) 557. 

[27] F. E. Close and R. G. Roberts, Phys. Lett. B336 (1994) 257. 

[28] See S. J. Brodsky, M. Burkardt and I. Schmidt, Nud.Phys. B441 (1995) 197. 

[29] R. L. Heimann, Nucl. Phys. B64 (1973) 429. 

[30] B.I. Ermolaev, S.I. Manayenkov and M.G. Ryskin, DESY preprint 95-017. 

[31] R. D. Ball and S. Forte, Phys. Lett. B335 (1994) 77. 

[32] S. Forte, Phys. Rev. D47 (1993) 1842. 

[33] B. Ehrnsperger and A. Schafer, Phys. Lett. B348 (1995) 619; 

J. Lichtenstadt and H. J. Lipkin, Phys. Lett. B353 (1995) 119. 

[34] B. Webber, summary talk at the International Conference on High Energy Physics 

(ICHEP), Glasgow, Scotland, 1994, [hep-ph/9410268. 



15 



parm. 


AB 


A T) 

An 


US 




0.48 ± 0.09 


0.42 ± 0.07 


0.35 ± 0.04 


% 


1.52 ± 0.74 


1.11 ± 0.52 


0.99 ± 0.23 


"NS 


—0.68 ± 0.15 


—0.67 ± 0.13 


—0.70 ± 0.15 


aq 


0.41 ± 0.38 


0.34 ± 0.59 


0.97 ± 0.65 


ag 


—0.47 ± 0.30 


0.03 ± 0.63 


0.35 ± 0.67 


Pns 


o o _i_ n o 
Z.Z ± U.O 


1 o _i_ n o 
l.O ± U.O 


1.6 ± U.O 


Pq 


O O 1 "1 /I 

3.3 ± 1.4 


2.1 ± 1.0 


r\ o \ "1 o 

0.8 ± 1.2 


Pg 


O 1 /I o 

Z.D ± 4.0 


( .U ± o.o 


4 (^nxeaj 




1 « _|_ 1 7 
ID It i / 


O.D It O.Z 


10 It iO 


aq = Gg 


0.1 ±3.0 


0.9 ±4.4 


-1.25 ±0.09 




62.6/63 


61.0/63 


61.8/64 



Table 1: Best-fit values of the parameters eq. ( p^ ) and for fits in the AB, AR and OS 
schemes; the errors shown are statistical only. 





Vq 


Vg 


r?(io) 


rf(io) 


ao(lO) 




as tab. 1 


0.48 ±0.09 


1.52 ±0.74 


0.122 ±0.013 


0.025 ±0.013 


0.14 ±0.10 


62.6 


high 6ns 

low 6ns 


0.48 ±0.06 
0.50 ±0.10 


1.43 ±0.48 
1.55 ±0.81 


0.124 ±0.011 
0.123 ±0.015 


0.027 ±0.011 
0.027 ±0.015 


0.15 ±0.08 
0.15 ±0.11 


63.2 
64.3 


ons = 


0.37 ±0.09 


1.17 ±0.70 


0.120 ±0.010 


0.023 ±0.010 


0.11 ±0.08 


73.9 


high gA 
low Qa 
high ag 
low as 


0.48 ±0.09 
0.48 ±0.09 
0.45 ±0.09 
0.51 ±0.09 


1.54 ±0.74 
1.50 ±0.73 
1.58 ±0.74 
1.47 ±0.73 


0.124 ±0.013 
0.120 ±0.012 
0.122 ±0.013 
0.122 ±0.012 


0.025 ±0.013 
0.026 ±0.012 
0.026 ±0.013 
0.025 ±0.012 


0.14 ±0.10 
0.14 ±0.10 
0.10 ±0.10 
0.18 ±0.10 


62.5 
62.7 
62.7 
62.5 


high as 
low as 


0.51 ±0.08 
0.47 ±0.09 


1.36 ±0.55 
1.80 ±0.82 


0.118 ±0.013 
0.124 ±0.011 


0.022 ±0.013 
0.026 ±0.011 


0.12 ±0.10 
0.14 ±0.09 


63.6 
62.1 


high fact, 
low fact. 


0.42 ±0.05 
0.41 ±0.07 


1.17 ±0.45 
0.81 ±0.39 


0.121 ±0.010 
0.133 ±0.007 


0.023 ±0.010 
0.037 ±0.007 


0.11 ±0.08 
0.25 ±0.06 


61.1 
63.0 


high ren. 
low ren. 


0.42 ±0.05 
0.42 ±0.08 


1.31 ±0.57 
0.92 ±0.59 


0.125 ±0.010 
0.129 ±0.016 


0.028 ±0.010 
0.033 ±0.016 


0.15 ±0.07 
0.20 ±0.13 


61.4 
62.3 


high thr. 
low thr. 


0.57 ±0.06 
0.46 ±0.07 


1.69 ±0.33 
1.54 ±0.67 


0.121 ±0.006 
0.121 ±0.012 


0.024 ±0.006 
0.025 ±0.012 


0.13 ±0.06 
0.13 ±0.10 


62.5 
62.0 



Table 2: Values of the parameters r/q and the first moment of gi eq. @, the axial charge 
eq. @ and (64 degrees of freedom for the entry ons = and 63 d.f. for all other en- 
tries), all calculated in AB scheme, for the various fits discussed in the text. 
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source 


Ar?(io) 


Arf(io) 


Aao(lO) 


fit (statistical) 


±0.013 


±0.013 


±0.10 


fitting procedure 


±0.003 


±0.003 


±0.03 


SU(2) violation 


±0.002 


±0.000 


±0.00 


SU(3) violation 


±0.000 


±0.000 


±0.04 


value of as 


+0.002 
-0.004 


+0.001 
-0.003 


+0.00 
-0.02 


thresholds 


±0.001 


±0.001 


±0.01 


higher order corrns. 


+0.011 
-0.002 


+0.012 
-0.003 


+0.11 
-0.03 



Table 3: Errors in the determination of Fi and ao- 



17 



1.50 



1.25 



1.00 



g?(x,Q') 0.75 



0.50 



0.25 



0.00 

10" 



"1 1 — I — r 



dashed: Q^ = Q^(x) (SMC) 
solid: Q^=Q^(x) (E143) 
dotdashed: Q^=l GeV^ 
dotted: Q^ = 10 GeV^ 



10" 




0.50 



0.25 



0.00 



g?(x,Q') -0.25 



-0.50 



-0.75 



-1.00 
10 



"1 1 1 — I — r-n 




■ / 



dashed: Q^ = Q^(x) (SMC) 
solid: Q^ = Q^(x) (E143) 



dotted: Q^=10 GeV^ 



-3 



10 



-2 



10 



-1 



lO*- 



Fig. li 



0.06 



0.04 



xg?(x.Q^) 



0.02 



0.00 

10 



dashed; Q^=Q^(x) (SMC) 
solid: Q^ = Q^(x) (E143) 
dotdashed: Q^=l GeV^ 
dotted: Q^=10 GeV^ 




10" 



10" 



10'" 



0.04 



xg?(x,Q') 0.02 



0.00 



10 



dashed: Q^ = Q^(x) (SMC) 
solid: Q^ = q^(x) (E143) 
dotdashed: Q^=l GeV^ 
dotted: Q^=10 GeV^ 




Fig. lii 

Plots of gi{x) (i) and xgi{x) (ii) compared to the SMC (crosses) and E143 (dia- 
monds) experimental data for (a) proton and (b) deuteron. The curves correspond to a 
NLO computation in the AB scheme with the initial parton distributions eq. ( p^) and 
the values of the parameters given in table 1. 



